Homogenization of the Vlasov equation and of the 
Vlasov - Poisson system with a strong external mag- 
netic field 



Emmanuel Frenod 

LMAM, Universite de Bretagne Sud 

Centre Universitaire de Vannes 

1 rue de la loi, F-56000 Vannes 

E-mail: Emmanuel.Frenod@univ-ubs.fr 

Eric Sonnendriicker 

CNRS, Institut Elie Cartan, Universite Henri Poincare - Nancy 1 
B.P. 239 

F-54506 Vandceuvre-les-Nancy Cedex 
E-mail: Eric.Sonnendrucker@iecn.u-nancy.fr 

Abstract: Motivated by the difficulty arising in the numerical simulation of the move- 
ment of charged particles in presence of a large external magnetic field, which adds an 
additional time scale and thus imposes to use a much smaller time step, we perform in 
this paper a homogenization of the Vlasov equation and the Vlasov-Poisson system which 
yield approximate equations describing the mean behavior of the particles. The convergence 
proof is based on the two scale convergence tools introduced by N'Guetseng and Allaire. 
We also consider the case where, in addition to the magnetic field, a large external electric 
field orthogonal to the magnetic field and of the same magnitude is applied. 

Keywords: Vlasov-Poisson equations, kinetic equations, homogenization, gyrokinetic 
approximation, multiple time scales, two scale convergence. 

Abbreviated title: Homogenization of the Vlasov equation 



f 



1 Introduction 



In many kind of devices involving charged particles, like electron guns, diodes or tokamaks, a 
large external magnetic field needs to be applied in order to keep the particles on the desired 
tracks. In Particle-In-Cell (PIC) simulations of such devices, this large external magnetic 
field obviously needs to be taken into account when pushing the particles. However, due 
to the magnitude of the concerned field this often adds a new time scale to the simulation 
and thus a stringent restriction on the time step. In order to get rid of this additional time 
scale, we would like to find approximate equations, where only the gross behavior implied 
by the external field would be retained and which could be used in a numerical simulation. 

The trajectory of a particle in a constant magnetic field B is a helicoid along the magnetic 
field lines with a radius proportional to the inverse of the magnitude of B. Hence, when this 
field becomes very large the particle gets trapped along the magnetic field lines. However 
due to the fast oscillations around the apparent trajectory, its apparent velocity is smaller 
than the actual one. This result has been known for some time as the "guiding center" 
approximation, and the link between the real and the apparent velocity is well known in 
terms of B. We refer to Lee [13] and Dubin et al [5 J for a complete physical viewpoint review 
on this subject. In the case of a cloud of particles, the movement of which is described by the 
Vlasov-Poisson equations, the situation is less clear as in the one particle case because of the 
non linearity of the problem. In particular, the question of knowing if the mutual influence 
of the particles can be expressed in terms of their apparent motion or if the oscillation 
generates additional effects is important. 

In this paper, we deduce the "guiding center" approximation in the framework of partial 
differential equations via an homogenization process on the linear Vlasov equation. Then, 
we show that in a cloud of particles the mutual influence of the particles can be expressed 
in term of their apparent motion without any additionnal terms. This is provided applying 
an homogenization process on the Vlasov-Poisson system, similar to the one used for the 
linear Vlasov equation and using the regularity of the charge density. 

Hence, we apply first a homogenization process to the linear Vlasov equation with a 
strong and constant external magnetic field. In other words for a constant vector J\A G i 2 , 
we consider the following equation: 



In this equation f £ = / e (t,x, v) with t G [0,T), for any T G M + , x G IR? X and v G 1R%. For 
convenience, we introduce the notations Q = M x x 2R^, O = [0, T) x M x and Q = [0, T) x Q. 
The initial data satisfies 



The electric field E £ = E e (t, x) and the magnetic field B e = B £ (i, x) are defined on O, both 
bounded in L°°(0, T; L 2 (Ml)) and satisfy 







E e — > E in L°°(0,T;Ll c (lRl)) strong, 
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and 

B £ — > B in L°°(0,T;LL(^)) strong, (1.4) 

for any T E M + . 

With those conditions, for every e there exists a unique solution of the equation (|l.ip 
/ e € L°°(0,T; L 2 (fi)). We characterize here the equation satisfied by the limit / (in some 
weak topologies) of (/ e ) e . 

The first main result of the paper is the following: 

Theorem 1.1 Under assumptions il.2\) - (L4\ ), the sequence (f £ ) E of solutions of the Vlasov 
equation satisfies, for any T £ M + , 

f e ->> f in L°°(0,T;L 2 (ft)) weak-*. 

Moreover, denoting for any vector v, Vii = (v • A4)A4, f is the unique solution of: 

Of 

+ V|| ■ V,/ + (Ey + v x B„) • V„/ = 0, 

dt 1 ,2* (1-5) 

f\t=o = 7T~ / /o( x > u ( v , T )) dr, 
2ir Jo 

where u(v, r) is i/ie rotation of angle t around A4 applied to v (see \2. 7|) /or more details). 



The deduction of this Theorem uses the notion of two scale convergence introduced by 
N'Guetseng [T5| and Allaire [2j. This convergence result is the following: 

Theorem 1.2 (N'Guetseng 115ij and Allaire If a sequence f e is bounded in L°°(0, T; i 2 (0)), 
then for every period 9 , there exists a 9 -periodic profile Fg(t,T,x,v) € L°°(0, T; L^(1R T ; L 2 (ft))) 
such that for all ipg(t,T,x,v) regular, with compact support with respect to (t,x, v) and 9- 
periodic with respect to r we have, up to a subsequence, 

( fipidtdxdv -»• / / FgTpgdTdtdxdv. (1.6) 

Above, Lf{IR T ) stands for the space of functions being L°°(M) and being 9-periodic and 

i>% = il>e(t, i,x,v). 

The profile F$ is called the 9-periodic two scale limit of f £ and the link between Fg and 
the weak—* limit f is given by 

/ Fg(t,T,X,v)dT = f(t,X,\). (1.7) 

Jo 

This result has been used with success in the context of homogenization of transport equa- 
tion with periodically oscillating coefficients by E [6] and Alexandre and Hamdache [lj and 
in the context of kinetic equations with strong and periodically oscillating coefficients in 
Frenod [7] and Frenod and Hamdache [8]. 

Here, since the strong magnetic field induces periodic oscillations in f £ , the two scale limit 
describes well its asymptotic behavior. Therefore, it is the right tool to tackle homogeniza- 
tion of equation (jl.lj) , and Theorem 11.11 is a consequence of the following result concerning 
the two scale limit of f £ . 
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Theorem 1.3 Under assumptions \l.ty) - [T^\ ), the 2n -periodic two scale limit F € L°°(0, T; 
L2^{M T \ L 2 (Cl))) of f £ is the unique solution of 

( dF 

— + (vxM)- V V F = 0, 

' — + v,| • V X F + (E|| +v x By) • V v F = 0, (1-8) 
dt 1 

- F]t=0 = ^^°( x ' u ( v ' r ^' 

where u(v,r) is the rotation of angle r around Ai applied to v (see {2.7$ for more details). 

Remark 1.1 This last Theorem can be generalised to the case of a non uniform strong 
magnetic field, see Remark 12.21 

Remark 1.2 The method we develop enables also to deduce the homogenized equations 
when a strong external electric field is added in equation (jl.ip . In this case we get, in 
addition aA^xM drift This drift shares a lot with the one used in the "guiding center" 
approximation. See Theorems 13.11 and 13.21 for more details. 

Secondly, as a relatively direct application of the former results, we may characterize the 
asymptotic behavior of a sequence (f £ ,W) of solutions of the following Vlasov-Poisson 
system 

( df £ M 

4- + v • V x f B + (E £ + v x — ) • V v f = 0, 

< f T f £ (1.9) 

J\t=o - Jo, y ' 

E e = -Vu £ , -Au £ = p £ , 

with 

P £ = ! J £ dv, 

where /o is assumed to satisfy 

/o > 0, /o G L 1 ni 2 (fi), and < / / (1 + \v\ 2 )dxdv < oo. (1.10) 

Jn 

The second main result of the paper is the following: 

Theorem 1.4 Under assumption $1.10\) . for each e, there exists a solution (/ e ,E e ) of M.9\) 
in L^O^L 1 n L 2 {tt)) x L°°(0,T; W 1 'i{M%) for any T G M + . Moreover this solution is 
bounded in L^^TjL 1 nL 2 (Sl)) x L°°(0,T; W l 'i(M^.)) independently on e. 

If we consider a sequence (/ e ,E £ ) of such solutions, extracting a subsequence, we have 

f £ -± f in L°°(0,T;L 2 (Q)) weak-*, 
E £ -> E in L°°(0,T;L 2 oc (iR3)) strong, 
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and, denoting for any vector v, vii = (v • M)M, the limit (/, E) satisfies 

r of 



+ V|, ■ V*/ + E|| •V«/ = J 
1 /" 27r 

/|t=o = ^ / /o(x, u(v, r)) dr, 
Vu, —An = p, 



94 
E 



(1.11) 



with 



P 



f fdv, (1.12) 
where u(v, r) is i/ie rotation of angle r around A4 applied to v (see 7| ) /or more details). 

This Theorem is a consequence of the following result. 



Theorem 1.5 Under assumption U.10\) . extracting a subsequence, the 2ir— periodic two 
scale limit F G L°°(0, T; L^(M T ; L 1 n L 2 (Q.))) of f £ is solution of 



( OF 

-^ + (vxM)-V v F = 0, 
dF 

— + V|, • V X F + E|, • V„F = 0, 
dt 1 

F\t=o = 2^/o(x,u(v,t)) 
E = — Vn, —An = p, 



(1.13) 



TOi/i /O given by 



Fdv. 



Moreover, p does not depend on t and 



p = p. 



(1.14) 
(1.15) 



This Theorem provides a rigorous justification of the procedures called subcycling and orbit 
averaging that are often used in Particle- In-Cell simulations of the Vlasov-Poisson equations 
in order to reduce the cost of the simulation. This procedure consists in advancing the 
particles, which provide an approximate solution of the Vlasov equation, on a smaller time 
step than the one used to advance the solution of the Poisson equations, see [3] for an 
overview. 



Remark 1.3 We may add in (jl.9p an external magnetic field B e strongly converging to B. 
The Theorems 11.31 and 11.11 will also apply in the same way leading to equation (|1 . 13[) with 
(|1.131 b) replaced by 

dF 

— + v,| • V X F + (E|, + v x By) • V„F = 0, 

and to (jl.lip with (jl.lli a) replaced by 

Of 



+ v,| • V x f + (E,| + v x B,|) • V v f = 0. 
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Remark 1.4 the results of this paper are connected to the ones of Grenier [91 [10] and 
Schochet [T7] concerning the pertubation of hyperbolic systems by a ^—depending linear 
operator inducing fast oscillations in time. In those work, the fast oscillations are canceled 
applying a "reverse-oscillating" operator to the familly of solutions of the e— depending 
equation whose limit is seeked. 

Yet, in the present work, the oscillations are treated using an ad-hoc class of oscillating test 
functions. 

From now on, and with no loss of generality, we set A4 = e±, where ei is the first vector 
of the basis (ei,e2,e3) of M^. 

The paper is organized as follow: In section [2j we provide the homogenization of the 
Vlasov equation (jl.ip . For this purpose, we first deduce the equation satisfied by the two 
scale limit F of f £ and then we obtain the homogenized equation of Theorem 11.11 
In section [3] we study the case when a strong electric field, orthogonal to the strong magnetic 
field, is added in the Vlasov equation. 

The last section is devoted to the homogenization of the Vlasov- Poisson system (jl.9j) . This 
is an application of the result proved in section 2] once the regularity of p £ implying the 
strong convergence of E £ is exhibited. 

2 Homogenization of the Vlasov equation with a strong ex- 
ternal magnetic field 

2.1 Two scale limit of the Vlasov equation 

Let us derive first the classical a priori estimate that are available for the Vlasov equation 

OH). 

Lemma 2.1 Under assumption h 1.2(1 on the initial condition /o and tl.3\ ), fcL4\ ) on the fields, 
there exists a constant C independent of e such that the solution f £ of the Vlasov equation 
hl.l\) satisfies: 

\\f £ \\L°°(0,T;L 2 (n)) < C, (2.1) 

Proof. We multiply the Vlasov equation (jl.ip by f £ and easily get 

4- I (/ £ ) 2 dxdv+ / wV^/fdxdvi / (E £ + vx (B £ + — )) • V4/ £ ) 2 dxdv = 0. 
at Jn Jn Jn e 

Integrating the second and third terms by parts, they vanish. Hence 

^ / (/ £ ) 2 dxdv = 0, 
at Jn 

which means that the 1? norm of f £ is conserved and gives us the estimate on f £ thanks 
to the bound on the initial condition /q. ■ 
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From the a priori estimate obtained in lemma [2TT1 we deduce that up to a subsequence, 
still denoted by e: 

f f in L°°(0,T;L 2 (ft)), weak-*. 

Applying the result of N'Guetseng [15] and Allaire [2] exposed in the Introduction (see 
Theorem II. 2p we may deduce that for every period 8, there exists a ^-periodic profile 
F e (t,T,x,v) G L oo (0,T;Lg o (iR T ;L 2 (^))) such that for all ^(t,r,x,v) regular, with com- 
pact support with respect to (t, x, v) and ^-periodic with respect to r we have, up to a 
subsequence, 



/ fiffidbdxdv-t I I Fgipgdrdtdxdv. (2.2) 

JO. J Q.Jo 



The two scale limit of equation (jl.ip is led in three steps. In the first one, using a 
weak formulation of (jl.lj) with oscillating test functions (see Tartar [18J and Bensoussan, 
Lions and Papanicolaou [3]), and using the convergence result (|2.2j) . we obtain a constraint 
equation for the ^-periodic profile Fg, for every 8. The second step is devoted to the 
resolution of this constraint equation. This will lead to the natural period 8 = 2tt for the 
profile. In the third one, using ad hoc oscillating test functions we deduce the equation 
satisfied by F. 

Step 1. Deduction of the constraint equation: Multiplying the Vlasov equation (jl.ip 
by ipg = tpg(t, ~,x, v) with ij)g(t,T, x, v) regular, with compact support in (i,x, v) G Q and 
^-periodic in r € M T , and integrating by parts using that div v (E £ + v x (B £ + =^)) = 0, we 
get 

/ rii^y + -(^r + v • (v^r + (e £ + v x (b £ + — )) ■ (v,v e ) £ ] a d^ = 

Jq at e or e 

- / fo(ipe)f t =o dxdv - 
Jn 1 

(2.3) 

Notice that Q = [0, T) x f2 is not an open set. Consequently, the contribution of the initial 
data stand in the left hand side of (|2.3p . 

Multiplying ()2.3[) by e and passing to the limit, applying (|2.2p . we deduce that the ^-periodic 
profile Fg associated with f £ satisfies 

dibg 

Fe(-^r + {vxM)-V v ^jg)dTdtdxdv = 0, ^ ^ 

Fg is 8— periodic in r. 

As this is realized for every ijjg regular, compactly supported in (t, x, v) G Q and ^-periodic 
in r G IR T , it is straightforward to show that (12. 4p is equivalent to 




/ / & Fg{^^- + (v x M.) ■ V v ipg) dr dt dv = 0, 

J Ml JO OT 



(2.5) 



Fg is 6*— periodic in r, 
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for every ipg = ^pe(r, v) regular, with compact support in v and ^-periodic in r, for almost 
every (x, t), and also equivalent to 

^ + (v x M) ■ V v F e = in V'(M T x iR^ 
is 9— periodic in r, 

for almost every (x, t). The two equivalent forms (|2,5p and (|2.6[) will be used in the sequel. 
Since the solution of this equation is not unique, we call it "constraint equation", and the 
goal is now to derive the form this equation imposes to Fg. 

Step 2. Consequences of the constraint equation: Intuitively, the constraint equation 
means that Fg is constant along the characteristics of the dynamical system V = V x Ai . As 
we shall see soon, those characteristics are helicoids around the magnetic vector Ai. Hence 
we first introduce u(v, r) the transformation v € M? — > u(v, r) € M? letting invariant the 
projection of v onto Ai and rotating of an angle r its projection onto the plane orthogonal 
to Ai. We have: 

u(v, r) = v\e\ + [v2 cos r — v 3 sin r]e2 + [v 2 sin t + v% cos T]e$. (2.7) 

Forgetting for the time being the periodicity ondition let us see what the constraint 
equation yields. Consider F(t,v) solution of the following equation 

dF 

— + (v x M) ■ V V F = in V'(M T x Ml). (2.8) 

Lemma 2.2 If F(r,v) G L°°(M T , L 2 (Ml)) satisfies 12.8]) , then there exists a function G € 
L 2 (M? U ) such that 

F(r,v) = G(u(v,r)). (2.9) 



Proof. Following Raviart [16], where weak solutions of first order hyperbolic equations 
are derived using their characteristic equations, F satisfies (|2.8p if and only if it is the 
translation along the characteristics V(r; v, s) solution of 

— = V x M, 

dr (2.10) 
V(s;v,s) = v, 

of a function G. This means 

F(r,v) =G(y(0;v,r)). 

Remark 2.1 In ()2.10j) V(t;v, s) means the solution at time t = r taking the value v at 
time t = s. 

As Ai is the first basis vector ei of the basis (ei, 62,63) of Ml an easy computation yields 
the solution of (j2~T0|) : 

V(r; v, s) = v±ei + [t>2 cos(t — s) + V3 sin(r — s)]e2 + [— V2 sin(r — s) + V3 cos(r — s)]e3. (2.11) 
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So we have 

K(0;v,r) = u(v,r), 

where u(v,r) is defined by (|2.7p . And F(r, v) = G(u(v,r)). ■ 

Now we take the periodicity condition under consideration. First, because of the previ- 
ous Lemma and since u(v, r) is 27r-periodic in r, if the period 9 is incommensurable with 
2-7T a fonction F(t, v) solution of (|2.8p and 0— periodic in r is constant. As a consequence, 
the 0— periodic profile Fg contains no information on the oscillations of (/ e ). 
Beside this, if 6 is (a multiple of) 2ir, a fonction F(t, v) satisfying (12.8f) naturally satisfy 
the periodicity condition. More precisely we have the 

Lemma 2.3 A function F(t,v) G L^(M t , L 2 (Ml)) satisfies 

dF 

— + (v x M) ■ V V F = in V'(M T x M? v ), 

if and only if there exists a function G G L 2 (Ml) such that 

F(r,v) = G(u(v,r)). (2.12) 

Hence, among every possible profile, we are incited to select the 2ir— periodic one. 

As a conclusion of this step, applying Lemma 12.31 the 2-7r-periodic profile F(t,r,x,v) = 
F2n(t,T,x,v) associated with the sequence f £ , solutions of (jl.ip . writes 

F(t,r,x,v) = G(i,x,u(v,r)), (2.13) 

for a function G G L°°(0,T; L 2 (1R 3 X x Ml)). 

Step 3. Equation satisfied by G: We now look for the equation satisfied by G linked to 
F by (I2TT3]1 . We have 

Lemma 2.4 The function G(t,x,u) linked to the lit -periodic profile F by &2.13]) is the 
unique solution of: 

dG 

— + u,| • V.G + (En + u x By) • V U G = 0, 

dt 11 1 11 11 (2.14) 

G,t=o = 2^/o, 

where for any vector field u = uiei + U2&2 + ^363, we define un = «iei its projection onto 
the direction of the external magnetic field A4. 

Proof. For any regular function tp, let us consider the regular and 27r-periodic in t function 
ip(t, r, x, v) = ip(t, x, u(v, t)) which satisfies 

^ + (vxM)-V v ip = 0. (2.15) 
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Now take ip e = tp(t, |,x,v) as a test function in the weak formulation (|2.3[) of the Vlasov 
equation (jl.ip . Then due to (|2.15p . we have 

J f £ [(%Y + v • (V^) £ + (E £ + v x B £ ) • (V^) £ ] dtd*dv = - f (^ t=0 dv. (2.16) 
Passing to the two scale limit, using that E e and B e converge strongly, yields 

/ / F{^-+vX7 x ip+(E+vxB)-(X7 v i))}dTdtdxdv = -[ f i/>(0, 0,x,v)dxdv, (2.17) 
JqJo at Jn 

which is equivalent to 

J G(t, x, u(v, r)) (t, x, u(v, r)) + v V x (p(t, x, u(v, r)) 

+(E + v x B) • V v {tp(t, x, u(r, v))}] dr dtdxdv = — [ / ¥>(0,x, u(v, 0)) dxdv. (2.18) 

Jn 

Making the change of variables u = u(v, r), noticing that v = V(r; u, 0) with the notations 
of (|2.10p and that du = dv, we get on the left-hand-side of (|2. 18jl 

f f 2w dip 

/ / G(i,x, u)[— (£,x, u) + [ttiei + [u2 cost + u 3 sin r]e 2 
JQJO OT 

+ [— ii2 sinr + U3 cos T]e 3 ] ■ V x tp(t, x, u) + (E + [u\ei + [112 cos r + u 3 sinr]e2 

r ■ 1 i n\ r r dtp . dtp . 

+ 1—U2 sinr + u 3 cosrJe 3 x B) • - — ei + [cost- h suit- — Je 2 

+ [-sinT-^- + cosT-^]e 3 ] drdtdxdu. (2.19) 

0U2 0113 

Now as neither G nor tp depend on r, let us perform the integration with respect to r. This 
yields, dividing by 2tt, 

/" c^c^ dip dip 

/ Gf— + ujei • VxV + (-Bit; 1" J3i(u 3 - u 2 t; — )] dtdxdu 

Jq at aui 0U2 ou 3 

= -— [ /ov?(0,x,u))(ixdu. 
27T Jn 

as the terms in cost, sinT and cost sin r vanish when integrated between and 2tt. This 
gives us the equation verified by G in the sense of distributions. ■ 

The uniqueness of the solution of (|2.14p enables to deduce that the whole sequence 
f £ two scale converges to F and, because of the link (jl.7p between F and /, weakly—* 
converges to /. 

2.2 Proof of Theorems (TTT] and Q 

Replacing first u by u(v, t) in ([2.14p we have 
dG 
~dt 



dG 

— (t,x,u(v,T)) + U||(v,t) • V x G(t,x,u(v,r)) 



+ (E||+u(v,r) xB||)-(V u G)(t,x,u(v,T)) =0, (2.20) 
G(0,x,u(v,t)) = ^-/ (x,u(v,t)). 
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Secondly, we get 

(E,| +u(v,r) x B||) • (V„G0(t,x,u(v,T)) = (Ey + v x By) • V«(G(t,x,u(v,T))). (2.21) 
Indeed, straightforward computations give 

Si — (G(i,x,u(v,r)))=S 1 — (t,x J u(v,r)), 

B1U3- — (G(t,x, u(v,r))) = -Bi(u 2 sin r cost- |-i>3Cos 2 r- h 

2 dG dG . 

f 2 sin T— hl>3 SUIT COST— — ), 

5 dG <9G 

B\U2t, — (G(t,x, u(v, r))) = — Bife sinrcosr- hu 3 sin 2 r- h 

0V2 0U2 0U2 

2 dG dG, 
V2 cos r— ^smrcosr— — J, 

and summing up these three relations yields (|2.2ip . Hence (|2.20p writes 
dG 

— (t,x,u(v,r)) + U||(v,r) • V x G(t, x, u(v, r)) 

+ (E[|+vxB||).V„(G(t,x,u(v,r)))=0 J (2.22) 
G(0,x,u(v,r)) = — / (x,u(v,r)), 

and since 

F(t,r,x,v) = G(t,x,u(v,r)), 
we obtain equation (|1.8l b.c). As (jl.8l a) is only the constraint, Theorem 11.31 is proved. ■ 

Then, in order to achieve the proof of Theorem[L3j we shall deduce the equation satisfied 
by / using the integral relation (|1.7p linking F and /. 

As 

/(t,x,v)= / F{t, t, x, v) dr = / G(t, x, u(v, r)) dr, 
Jo Jo 

integrating (12.22P in r, gives 

^(J G{t,x t u(y,r))dTj+v ]r V x (J G(i,x, u(v,r)) dr 

(E|| + v x B||) • V„ ^ G{t, x, u(v, r)) dr^j = 0, (2.23) 

/•'An Y f^K 

/ G(0,x,u(v,r))dr = — / / (x, u(v, r)) dr, 
Jo 27T Jo 



G(0,x,u(v,r))dr = — / 

/0 27T Jo 

giving the homogenized Vlasov equation (jl.5p . Hence Theorem 11.11 is proved 



Remark 2.2 The goal of this remark is to show how we may adapt our method to deduce 
the equation for the two scale limit in the case of a non uniform strong magnetic field. We 
assume that Ad(t, x) is a smooth function from O to the set of vectors with norme 1. That 
implies 

ex=K(t,x)M(t,x), (2.24) 
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where TZ(t, x) is a smooth map from O to the set of orthogonal matrices. In this framework, 
the way to deduce the constraint equation 



dF 

— + (vxM)-V v F = 0, 
or 



(2.25) 



remains the same. As 7£(i,x) is nothing but the matrix of a change of coordinates leading 
A4(t,x) onto ei this constraint implies 



F(t, t, x, v) = G(t, x, K T (t, x)u(Tl(t, x)v, r)) 



(2.26) 



where 1Z T is the transposed (but also the reverse) matrix of 1Z. In order to give a more 
usable form to (|2.26j) , we denote by 



r(r) 



/ 1 

cos r — sin r 
y sinr cost 



and since TZ T u(lZv, r)) = 1Z T r(T)7lv, (|2.26p also writes 

F(t,T,x,v) = G(t,x,^ T (t,x)r(r)^(t,x)v). 



(2.27) 



(2.28) 



Now, in order to get the equation satisfied by G, we proceed in a similar way as in the proof 
of Lemma I2.4I with this difference that we use test functions 



ip(t, t, x, v) = (p(t, x, TZ T (t, x)r(r)7£(i, x)v). 



(2.29) 



Those test functions naturally satisfy the constraint and the computation leading to equa- 
tion (I2.17P remains valid. But here, 



dip dip .dU T 



dt 

dip dip 



<)t + ^" d! 



dlZ 

r(r)ftv + n T r( T )— v) • V uV 
dt 



dxi 



,dll T , , m rp dlZ- 

_J_ / — i — ■ i ' ■ ! * I 

dxi dxi 



r(r)ftv + K T r{T)— v) • V u tp, 



(2.30) 
(2.31) 



and 



V v ip = {1l T r{T)K) T V u ip = {K T r(-T)K)V u ip. (2.32) 

Hence, in place of formula (|2.19p . making the change of variables u = 1Z t t{t)1Z\ , (v = 
1Z T r(— t)1Z\i), we get 



J jT G(t, x, u)[^ + (^TZu + n T r(r)^TZ T r(-r)nu) ■ V u 
+K T r(-T)Ku ■ (V x ip + 



( (^nn + TZ T r(T)^-TZ T r(-T)TZu)-V u p\ 

T 



+ n T v{r)^n T v{-r)nn) ■ V u p 



) 



\ (^u + ^r(r)^r(-r)7eu).V^ / 
+E • TZ T r(-T)TZV u ip + ({K t y(-t)1I\i) x B) • (K t y(-t)11V u cp)} dr dt dx du 

= - [ fo<p(0,x,u))dxdu. (2.33) 
Jo, 
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Since 

and 

with 



E • TZ T r(-T)TZV u ip = n T r{T)lZE ■ V u ip, 
{{1l T r(-T)Ku) x B) • (ft T r(-r)ftV u ^) = 1l T r(T)'llB'R T r(-T)Ku ■ V u <p, 



B 



B 3 -B 2 
-B 3 Bi 
B 2 -Bi 



we finally get, integrating (|2.33p in r and defining 

/ 1 








dG 

at 



+ TZ T rJZu ■ V X G + {K T TKE + — / * K T r{T)1lB'R T v{-T)'RdT u) • V U G 

2vr Jo 

+V U ■ [(^Ku + jf n T r(r)^n T r(-r)ndr u)G] 



+V M .[(f K T v{-r)nvi- 

27T Jo 



( (&u + ^r(r)g^r(-r)7lu)G^ 



V (ff7iu + ^r(r)g^r(-T)^u)G J 



G lt=0 - —fo- 



)dr] = 0, 

(2.34) 



where for 4 vectors of M 3 , A,B,C,D we denote 



A ■ 



( B 

C 



A l B + A 2 C + A 3 D. 



This last equation is the equation for the two scale limit when the strong magnetic field is 
non uniform. 



3 Homogenization of the Vlasov equation with strong exter- 
nal magnetic and electric fields 

We study here a variant of the previous problem. We homogenize the following Vlasov 
equation with strong and constant external magnetic and electric fields: 

df e M M 

+ v • V x f + ((E £ + T ) + v x (B £ + '—)) • V v f = 0, (31) 

f\t=o = /o> 
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for constant vectors M.€^ 2 ,Af^^ 2 ,A4-LJ\f and under the same assumptions (|1.2[> - (|1.4|) 
as previously. With no lost of generality, we set M. = e\ and M = e2- 

Since the a priori estimate (|2.ip remains valid we always have 

f -± f in L°°(0,T;L 2 (O)), weak-*. 

And, there exists a 27r-periodic profile F(t,r,x,v) G L°°(0, T; L^(M T ; L 2 (tt))) such that 
for all r, x, v) regular, with compact support with respect to (i,x, v) and 27r-periodic 
with respect to r we have 

f2w 



f f £ 4> 6 dtd-xdv -> f I Fipdrdtdyid-v. (3.2) 
JO. J Q.Jo 



IQ JQJO 

Proceeding as in section 12.11 we obtain the following weak formulation with oscillating 
test functions: 



r[( ^ + I ( ^ + v . (Vl „ 

7V\ _ M 



+ ((E £ + — )+vx (B £ + — ))-(V„V) e ]^x(iv = - / fo(ip)f t=0 dxdv, 
e e Jn 1 



(3.3) 



which, passing to the limit in e yields 



— + (M + v x M) ■ V V F = in V'(M T x Ml). (3.4) 



This constraint equation means that 

F(t,r,x,v) =G(t,x,y(0;v,r)) 
where V[t; v, s) is solution of 

dV 



M + VxM, , 
dr (3.5) 

F(s;v,s) = v. 



We have 



V(t; v, s) = viei+[t; 2 cos(r-s)+(v 3 +l) sin(r— s)]e 2 +[-U2 sin(r— s)+(v 3 +l) cos(r— s)— l]e 3 , 

(3.6) 

and then 

V(0; v, r) = u(v, r) = v \e\ + [v2 cos r — (^3 + 1) sin r]e2 + [+^2 sin r + (i>3 + 1) cos r — l]e3. 

(3.7) 

Hence we conclude 

F(t,r,x,v) = G(t,x,u(v,r)), 
for a function G G L°°(0,T; L 2 (lRl x 
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Remark 3.1 If M were not orthogonal to yVf, i.e. if N = n\e\ + e2, the consequence 
of equation (|3.4p would be F(t,r, x, v) = G(t,x, u(v, r) + nirei). Then the periodicity 
condition and the L 2 in v regularity of F would imply G = F = 0. 



Now in order to deduce the equation G satisfies, we take any regular function ip and we 
consider the regular and 2-7r-periodic in r function tp(t, r, x, v) = (p(t, x, u(v, r)) which 
satisfies 

^ + (^ + vxAf).V^ = 0. (3.8) 

Hence using this test function in the weak formulation (|3.3p . the term containing the con- 
straint disappears and passing then to the limit gives 

/ F[^- + v-V x 'ip + (E+vxB)-(V v i[>)}dTdtdyLdv = - [ f xp(0, 0, x, v) dxdv. (3.9) 
JqJo ot Jn 

Replacing F and ip by theire expressions, (j3.9j) becomes 

/q i) 2 " X ' U(V ' T)) ^ (t ' X ' U(V ' T)) + V ' Vx(p & x, u(v, r)) 

+(E + v x B) • V v {(p(t,x,u(T,v))}]dT dtdxdv = - [ f <p(0,x, u(v, 0)) dxdv. (3.10) 

Jn 

Making the change of variables u = u(v,r), we get on the left-hand-side of (|3.10|) 

f f 27T d(p 

/ / G(t,x, u)[— (t,x, u) + [u^ei + [u 2 cost + (u 3 + 1) sin r]e 2 

JQJO 

+ [-ii 2 shiT + (u 3 + 1) cost - l]e 3 ] ■ V x (p(t,x, u) 

+(E + [uiei + [u 2 cos t + (u 3 + 1) sin r]e 2 

+ [—M2 sin t + (tt 3 + 1) cos t - l]e 3 l x B) • [- — ei + [cos r- h sin r- — ]e 2 

dip dp 

+ [— sin t— h cos t— — ]e 3 ] dr dt dx du. 

OU2 ou 3 

(3.11) 

Now as neither G nor tp depend on r, let us perform the integration with respect to r. This 
yields, dividing by 2ir, 

f G[§ + (nie: - e 3 ) • + ((£i - S 2 )|^ + £i((«3 + 1)|^ - u 2 |^)] dxdu 

JQ Ot OU\ OU2 OU3 

/o¥>(0, x, u)) dxeiu. 



1 

"torJn 



Hence we proved the following result 



Theorem 3.1 Under assumptions U.^) - [L4^ and with Ai = ei and M = e 2 , the function 
G{t, x, u) linked to the 2ti -periodic profile F associated with the solution of (3. 1\) is the 
unique solution of: 



ui 



* \ 

G\ t=0 = —fo, 



V X G + 



1 





• V U G = 0, 



(3.12) 
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Now we deduce the equation satisfied by the weak limit / using the relation linking F 
and /. 

Replacing u by u(v,r) in A3. 12H . we have 



dG 
~dt 



m(v, t) 



(t,x,u(v,r))+ 



+ 

G(0,x,u(v,r)) 
Then as 




V x G(t,x,u(v,r)) 



/ ni(v,r) 
u 2 (v,t) 
V «3(v,r) + 1 




2tt 



/ (x,u(v,r)). 



• V u G(t,x,u(v,r)) =0, 



(3.13) 



V„G(t,x,u(v,r)) 




and as 



we have the 



V B (G(t,x,u(v,r))), 

(3.14) 

/(i,x,v) = F(t,r,x,v)dr= G(t,x,u(v,r))dr, 
Jo Jo 



Theorem 3.2 Under assumptions {H.2\ )- (L4\ ) and with A4 = e\ and N = e2, the weak—* 
limit f of the sequence of solutions of \3. 1\) is solution of: 




■ V*/ + 

1 f 2n 

f\t=o = 2^ J Q /o(x, u(v, t)) dr, 
with u(v, r) given by (3? 
Remark 3.2 As 




vi 

V2 



! 



\v 3 + 1 j ^ 



v,/ = 0, 



(3.15) 



MxM, 



-1 



the presence of —1 on the third componant of the advection vector shares a lot with the 
additional Af x J\A drift effect usually found in the "guiding center" approximation with 
same order electric and magnetic fields. 
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4 Homogenization of the Vlasov-Poisson system with a strong 
external magnetic field 



Assume that the initial condition satisfies the hypothesis (jl.lOp . Then for e > given there 
exists a solution in L°°(0, T; L 1 n L 2 (n)) x L°°(0,T; W l, i(IP? x )) for any T € M + of the 
Vlasov-Poisson system (ll.9p . Indeed for a given e > a force field associated to a constant 
magnetic field has simply been added to the system studied by Horst and Hunze [12] and 
their result can be extended easily to our case. 

Now in order to apply the results of section [2] and pass to the two-scale limit of the 
Vlasov equation, we need to prove that there exists a sequence of electric fields solution of 
the Vlasov-Poisson system (|1 .9H verifying 

E £ E in L°°(0,T;LL(^)) strong, 
for any T E M + . For this, we shall need some apriori estimates independent on e. 

Lemma 4.1 Assume that the initial condition /o is such that /o € L 2 (Q) and Jq fo\v\ 2 dxdv 
is bounded, then there exists a constant C independent of e such that the solution (E £ ,f £ ) 
of the Vlasov-Poisson equations hl.9\) satisfies, for any T € IR + , 

\\f £ \\L°°(0,T;L 2 (Q)) < C, 
\\(\v\ 2 n\\L°°(0,T;mn))<C, 

and moreover 

||p £ (x,i)|| 7 ro3 <C, 

||J £ (x,t)|| 7 ro3 < C, 

where p £ (x,t) = f f e dv and J £ (x, t) = f vf £ dv. 

Proof. Let us proceed formally in order to simplify the presentation. All the following is 
rigorously verified for regularized solutions and then the bounds, that do not depend on e, 
are conserved when passing to the limit with respect to the regularizing parameter. 

Multiplying the Vlasov equation by |v| 2 and integrating with respect to x and v, we get 
d 



/ / £ |v| 2 dvdx-2 / J e -E £ dx = 0. (4.1) 



dt 

Now, integrating the Vlasov equation with respect to v yields the continuity equation 

dp- 



m + V • J £ = 0. (4.2) 



Using this, we obtain 



J £ E £ dx = - J £ -Vu £ dx= V • J £ u £ dx, 
ml Jul Jnl 

dpe £ - 

-—- u dx. 

Ml dt 
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But the Poisson equation yields 

II / (Vn £ ) 2 dx= / ^u £ dx. 
2dtJml J Ml dt 

Hence (|4.ip becomes 

4(/ / e |v| 2 (ivdx+ / (Vu £ ) 2 dx) = 0. 
ere jo Jml 

Finally, integrating this last equation with respect to t and using the hypothesis on the 
initial condition /o, we get the bound on Jq / £ |v| 2 <ivcix. 

Let us now proceed with the L 2 estimate on f £ . To this aim we multiply the Vlasov 
equation by f £ and easily get 

4- I (/ E ) 2 dxdv+ / v V x (f £ ) 2 dxdv + / (E e + vx — )-V v (/ £ ) 2 dxdv = 0. 
at jo jo Jo e 

Integrating the second and third terms by parts, they vanish. Hence 

|/ n (/.)^v = 0, 

which means that the L 2 norm of f £ is conserved and gives us the estimate on f £ thanks 
to the bound on the initial condition Jq. 

We now come to the last estimates, namely those on p £ and J £ . Following the idea of 
Horst [11], the estimates for p £ and J £ can be obtained by decomposition of the velocity 
integral: 




for any R > 0. Using the Cauchy-Schwartz inequality, we have 

/ fdv<([ (/ £ ) 2 dv)l(/ dv)3 <di2f(/ (f £ ) 2 dv) l i, 

J\v\<R J\v\<R J\v\<R JMl 

and 

/ fdv<(f ^fdv<l/ M 2 fdv. 

J\v\>R J\v\>R R z R z Jml 

Hence, we have for any R > 

\pF(x,t)\ < dR$( [ (ffdvfi + -L f \v\ 2 fdv. 

J iRy ^ >J iRy 

Taking the R which minimizes the right-hand-side we obtain 

|p £ (x,t)| < C 2 ([ . (/ £ ) 2 dv)f ( / \v\ 2 f £ d^ , 
Jm 3 v Jml 
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and finally 

/ |p £ (x,t)|idx < CJ (/ 3 (/ £ ) 2 dv)§(/ H 2 / E dv)ldx, 

< c 3 (/ a (/ e ) 2 dxdv)f(/ M 2 rd X( fv)t, 

Jmlxml Jmlxml 

thanks to the Holder inequality. Now, knowing that the terms on the right-hand-side are 
bounded, we have our estimate on p e . 

We can proceed with J e in exactly the same way: 

|J £ (x,t)|< [ \v\fdv= [ |v|/ e dv + / |v|/ £ dv, 
Jml J\v\<R J\v\>R 

for any R > 0. Using the Cauchy-Schwartz inequality, we have 

/ |v|/ £ dv<(/ (/ £ ) 2 dv)!(/ j? 2 dv)i, 

J|v|<R J|v|<i? J\v\<R 

and 
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/ |v|/ e dv< / ^-fdv. 
J\v\>R J\v\>R it 



' £ + V • J e = 0. 



Here again we can find the minimizing R and use Holder's inequality to obtain 

/ |J £ (x,i)|£dx < C 4 ( / (/ e ) 2 dxdv)s(/ H 2 / £ dxdv)l, 
jjr3 J mix ml J mix ml 

which gives us the estimate on J E . ■ 

Now, on the one hand, thanks to the classical regularizing properties of the Laplacian, 
p e (x, t) bounded in L°°(0, T; Ls (IR^.)) implies that u £ such that — Au £ = p £ is bounded 
in L°°{0 } T-W 2 'l{lRl)) and hence E £ = -Vu e is bounded in L°°(0,T; W 1 '* . On 
the other hand, integrating the Vlasov equation with respect to v, we get the continuity 
equation 

d£_ 

dt 

Hence, as J £ is bounded in L°°(0, T; lI(MD), is bounded in L°°(0,T; W -1 '^ (M^.)), and 
as we have 

chx e _ <9p E 
~dT = ~dt' 

the regularizing properties of the Laplacian now yield ^ bounded in L°°(0,T; VF 1 '6(iR^)), 
which yields that 2|L is bounded in L°°(0, T; lI (JR-l)). 
Thus as 

W^ilRl) C Ll c (Ml) C LfjMl), 

the first injection being compact and the second being continuous, the Aubin-Lions Lemma 
(see for example Lions [2]) yields that the functional space 
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provided with the usual product norm, is compactly embedded in L°°(0,T; Lf oc (M^.)) and 
consequently, as E £ is bounded in U there exists a subsequence of E £ which converges 
strongly in L°°(0, T; Lf oc (Ml)). 
Hence we have 

Lemma 4.2 Under assumption ^1.10\) . extracting a subsequence, the sequence (f £ ,E £ ) G 
L^iO^-L 1 nL 2 (Q)) x L°°(0,T; W l, i{]R^)) of solutions of {Op satisfies: 

f -± f in L°°{0,T;L 2 (n)) weak-*, 

f two scale converges to F G L°°(0, T; L^(M T ; L 1 D L 2 (n))), 
E e -> E in L°°(0,T;Ll c (lRl)) strong. 

Moreover, 

p e p in L°°(0,T;Li(lRl)) weak-*, 
p e two scale converges to p G L°°(0, T; L^{M T ; Li (M^.))), 

J £ J m L°°(0,T;Li(iR3))) weafc- *, 
J e too sca/e converges to J G L°°(0, T; L^{M T ; (IR^))), 
/or any T G iR + . 

Then, passing to the two scale limit in (|1.9p applying the results of section [5] yields system 

In order to show that ^ does not depend on r, we multiply the continuity equation (14.2j) by 
tp(t, \,x) with </?(i, r, x) regular, with compact support in t, x and periodic in r, we integrate 
it by part and we multiply it by s. Passing then to the limit gives 




Now, summing up the above arguments proves Theorem 11.51 ■ 
Lastly, exactly as in section [21 by an integration with respect to r we get Theorem 11.41 ■ 
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